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Abstract 

We report measurements of partial branching fractions of inclusive charmless semilep- 
tonic 5-meson decays at the endpoint of the electron momentum spectrum. The 
measurements are made in six overlapping momentum intervals that have lower 
bounds ranging from 1.9 GeV/c to 2.4 GeV/c and a common upper bound of 
2.6 GeV/c, as measured in the centre of mass frame. The results are based on 
a sample of 29 million BB pairs, accumulated by the Belle detector at the KEKB 
asymmetric e + e~ collider operating on the T(4S I ) resonance. In the momentum 
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interval ranging from 1.9 GeV/c to 2.6 GeV/c we measure the partial branch- 
ing fraction AB(B -> X u ev e ) = (8.47 ± 0.37 ± 1.53) x 1CT 4 , where the first er- 
ror is statistical and the second is systematic. A prediction of the partial rate 
R = (21.69 ± 3.62^ i'gg) iK^pps -1 in this momentum interval based on theory is 
calculated with input HQET parameters that have been derived from Belle's mea- 
surement of the B — > X s j photon energy spectrum, where the first error is due to 
the uncertainty on HQET parameters and the second error is from theory. Using 
both AB{B -► X u eu e ) and R we find \V ub \ = (5.08 ± 0.47 ± 0.42 ±g;|g) x 10" 3 , where 
the first error is from the partial branching fraction, and the second and third errors 
are from uncertainties in R. 

Key words: Semileptonic i?-meson decays; CKM matrix 
PACS: 11.30.Er, 13.20.He, 12.15.Ff, 14.40.Nd 



1 Introduction 

The magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\Vub\ is a fundamental parameter of the Standard Model (SM). A knowledge of 
its value is crucial to the understanding of CP violation within the SM, which 
is underpinned by knowledge of the so-called Unitarity Triangle (UT). Recent 
precise measurements of UT angle <f>i(/3) [1,2] have focussed attention on \V u b\, 
because it determines the side of the UT that is opposite <f>i, it directly tests 
the Kobayashi-Maskawa mechanism [3] for CP violation within the SM. 

To date, inclusive measurements of \V u b\ have been reported by experiments 
operating on the T(4S) resonance, namely CLEO [4,5,6], ARGUS [7], BaBar [8] 
and Belle [9] , and by LEP experiments operating on the Z resonance, namely 
L3 [10], ALEPH [11], DELPHI [12] and OPAL [13]. 

The value of \V u b\ can be extracted from the measured rate of charmless 
semileptonic 5-meson decays B — > X u ev e in a kinematic region that has to 
be chosen to minimise the impact of the large background from the charmed 
semileptonic 5-meson decays B — > X c eu e . One such region is at the endpoint 
of the lepton momentum spectrum: in the rest frame of the decaying B meson, 
leptons from B — > X c ev e decays attain a maximum momentum of 2.31 GeV/c 
while for B — > X u eu e decays the maximum is 2.64 GeV/c. 

In this paper we report measurements of partial branching fractions to inclu- 
sive charmless semileptonic 5-meson decays from the yield of electrons and 
positrons in six overlapping momentum intervals. The intervals have lower 

1 on leave from Nova Gorica Polytechnic, Nova Gorica, Slovenia 



4 



limits commencing at 1.9GeV/c through 2.4GeV/c incremented in steps of 
0.1 GeV/c and a common upper limit of 2.6 GeV/c, as measured in the centre 
of mass (CM) 2 frame. 

We use two methods to extract \V ub \, one that has been standard practice [14] 
(DFN), and one that has been recently developed [15,16,17, 18, 19] (BLNP). 
The DFN method involves extrapolation from a partial to a full branching 
fraction using the DeFazio-Neubert prescription with given shape function 
parameters [20] followed by a routine to translate the full branching fraction 
to \V ub \ [14]. The BLNP method in contrast to the DFN method provides 
a more systematic treatment of shape function effects by incorporating all 
known contributions, includes power corrections, uses an improved perturba- 
tive treatment and directly relates the B — > X u eu e partial rate to F° r 
both methods we use values for the parameters of the shape function that 
were determined using the B — > X s j photon energy spectrum measured by 
Belle [21] 3 . In the case of DFN the procedure and measurements are given 
in Ref. [24]. The same procedure was updated with predicted shapes of the 
B — > X s 7 photon energy distributions from Ref [25] to yield values of shape 
function parameters relevant to the BLNP method. The latter are equivalent 
to HQET parameters in the shape function scheme [19]. 



2 Detector and Data Sample 

The results reported here are based on data collected with the Belle detector at 
the KEKB asymmetric energy e + e~ collider [26]. The Belle detector is a large- 
solid-angle magnetic spectrometer that consists of a three-layer silicon vertex 
detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- flight 
scintillation counters (TOF), and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect K\ mesons and to identify muons (KLM). The 
detector is described in detail elsewhere [27]. We use 27.0 fb _1 and 8.8 fb _1 
integrated luminosity samples taken at (ON) and 60 MeV below (OFF) the 
T(45 l ) resonance energy, respectively. The ON sample consists of 29.4 million 
BB events. Unless explicitly stated otherwise, all variables are calculated in 
the CM frame. 



2 The rest frame of the T(4S) is equivalent to the centre of mass frame. 

3 The use of the photon energy spectrum from inclusive radiative S-meson decays in 
determining the 6-quark shape function or distribution function was first discussed 
by Bigi et al [22] and Neubert [23]. 
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3 Data Analysis 



The procedure for this analysis largely follows that of CLEO [6] , and consists 
of examining the spectrum of electron candidates with momentum in the range 

1.5 — 3.5 GeV/c, which includes both signal and sideband regions. We initially 
chose and optimised our selection criteria for the momentum region, 2.2 — 

2.6 GeV/c. For ease of explanation we discuss the experimental procedure for 
this momentum interval and later describe the slight differences in the signal 
extraction for other momentum intervals. 

In the CM frame the kinematic endpoint for decays of the type B — > X c eu e , 
including the non-zero B momentum and detector effects, is 2.4 GeV/c. The 
B — > X u ev e signal is extracted in the momentum region 2.2 — 2.6 GeV/c (HI), 
while a lower range, from 1.5 — 2.2 GeV/c (LO), is examined to evaluate the 
contribution from B — > X c eu, which is then extrapolated to the HI region. 

The uncertainty on the fraction of B — > X u ev within the HI region is a ma- 
jor source of systematic error for the determination of the branching fraction 
and \V u i\. For choosing and optimising selection criteria we use a sample of 
events containing B — > X u ev e decays, generated via Monte Carlo simula- 
tion and based on a model described in Ref. [28] (ISGW2), which predicts 
the form factors and branching fractions of the many exclusive charmless 
semileptonic 5-meson decay channels that form the sample. We also gen- 
erate samples based on an inclusive B — > X u ev e model, according to the pre- 
scription of DeFazio and Neubert [20], with shape function (SF) parameters 
that correspond to the residual B- meson mass 4 and the average momentum 
squared of the 6-quark inside the 5-meson set to A SF = 0.659 GeV/c 2 and 
—\f F = 0.400 GeV 2 /c 2 , respectively. These values were determined from the 
photon energy spectrum in B — > X s 7 decays measured by Belle [24]. To exam- 
ine the extent to which our results may vary due to uncertainties in A SF and 
A SF , we also generate four samples with parameters that define the long and 
short axes of the A% 2 = 1 contour in the ( G ^v/ e 2 , gcvVc 2 ) P^ ane ' corresponding 
to (0.614,-0.231), (0.736,-0.714), (0.719,-0.462) and (0.635,-0.483). 

3.1 Event Selection 

To select hadronic events we require the multiplicity of tracks to be greater 
than four and the primary event vertex to be within 1.5 cm radially and 3.5 cm 
longitudinally from the detector origin. We make further requirements based 

4 A SF = M B - m b ; 

where rrif, and Mb are the masses of the 6-quark and i?-meson, respectively. 
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on quantities calculated in the CM frame - that the sum of cluster energies 
in the ECL satisfies 0.18\/s < -Eecl < 0.80^ where ^/s is the CM collision 
energy, that the visible energy E vis be at least 0.50^, that the absolute sum 
of longitudinal track and photon momenta be less than 0.30a/s, that the heavy 
jet mass be either at least 0.25 x E vis or greater than 1.8 GeV/c 2 , and that 
the average cluster energy be less than 1 GeV. We also require that the ratio 
i?2 of the second to the zeroth Fox- Wolfram moment [29] be less than 0.5. 



3.2 Electron spectrum 

Charged tracks are reconstructed from hits in the SVD and the CDC. Tracks 
are required to pass within a short distance from the interaction point (IP) 
of the e + e~ collision, where the £>-meson decays promptly. For improved data 
and MC agreement, tracks must be within the acceptance of the barrel part 
of the ECL; —0.63 < cos#i a b < 0.85, where #i a b is the polar angle measured 
in the laboratory frame with respect to the direction opposite to that of the 
positron beam. Tracks are identified as electrons on the basis of a matching en- 
ergy cluster in the ECL, and, subsequently, upon the ratio of ECL-measured 
energy to CDC-measured track momentum, transverse shower shape in the 
ECL, ionisation energy loss in the CDC, and the ACC light yield [30]. Given 
the track requirements, electrons with momenta in the range 1.5 — 2.6 GeV/c 
are positively identified with a probability of (94.0 ± 1.5)% while pions are 
misidentified as electrons with a probability of (0.13 ± 0.01)%, as measured 
using samples of reconstructed J/ip — > e + e~ and Kg — > it + it~ decays, respec- 
tively. 

To reduce the contribution of electrons from J/ip and ip{2S) decays and pho- 
ton conversions, our candidate electrons are paired with oppositely charged 
tracks identified as electrons in the event and rejected if their mass falls 
within either J/ip, ip(2S) or 7 mass windows, defined as [3.025, 3.125] GeV/c 2 , 
[3.614,3.714] GeV/c 2 and [0,0.1] GeV/c 2 , respectively. The photon conversion 
veto has the additional effect of removing electrons from tt° Dalitz decays. 
The yields of candidates that do not pass the J ftp veto requirement are com- 
pared in data and MC to determine a normalisation factor for MC-estimated 
backgrounds, which are described below. 

Since the dynamics of the hadronic part of B — > X u eu e is not well estab- 
lished, it is important that selection requirements retain acceptance over a 
wide range of q 2 = (p e + p u ) 2 (dilepton invariant mass squared) in order to 
minimise model dependence. Additional event requirements are designed to 
reduce contributions from continuum (e + e~ — > qq where q — u,d, s, c) and 
QED-related processes (including two photon and tau-pair events) without 
introducing a q 2 bias. A set of "energy flow" variables is formed by grouping 
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detected particles in bins of 0.05 in cos#, where 9 is the particle angle with 
respect to the candidate electron, and taking the energy sum in each bin. The 
energy flow in the backward direction —1.00 < cos# < —0.95 is not used, as 
it is found to disproportionately reduce the acceptance in the low q 2 region. A 
Fisher discriminant, denoted J^flow, is constructed from the remaining energy 
flow variables with coefficients chosen to best separate signal from continuum 
events. We also make use of a 6-quark rare decay tag variable, denoted /C, and 
calculated as 



where Q(e) is the charge of the candidate electron, and iV(i^ ± ) are the number 
of tracks identified as positively and negatively charged kaons in an event, 
respectively. K, exploits the presence of lepton-kaon charge correlations evident 
in BB events wherein one 5-meson decays via a b — > uev e transition whilst the 
other £>-meson decays typically via b — > c — > s transitions, thereby resulting 
in, on average, a net strangeness or kaon charge that is correlated to the charge 
of our candidate electron. The correlation does not exist in continuum events 
nor in BB events that do not involve the b — > uev e transition. Charged tracks 
are identified as kaons by utilising specific ionisation energy loss measurements 
made with the CDC, light yield readings from the ACC, and time of flight 
information from the TOF. The average kaon identification efficiency and fake 
rate in the momentum range 0.5-4.0 GeV/c, as measured in the laboratory 
frame, are (88.0 ± 0.1)% and (8.5 ± 0.1)%, respectively. 

To preserve signal efficiency, the selection requirements on JF FLOW are chosen 
differently for three cases of /C: /C > 0; /C = 0; and /C < 0. The cut values 
are chosen to optimise the figure of merit Sj \/S + B, where S and B reflect 
the signal and continuum background expectation, respectively, as estimated 
from MC events, assuming the branching fraction measured by CLEO [6] . The 
/C dependent cuts on .Fflow reduce continuum backgrounds by 97% while 
retaining 33% of the B — > X u eu e signal. 

To further suppress QED-related continuum backgrounds, the cosine of the 
angle between the thrust axis of the event and the e~ direction cos$£/j r , is 
required to be less than 0.75. Crucially, the thrust axis calculation includes the 
missing momentum as a component. Missing momentum is calculated as the 
difference between the momentum of the beams and the sum of the observed 
track and cluster momenta. Placing a constraint on cos 9 thr was found to bias 
the q 2 distribution in signal events less than a constraint imposed on the 
direction of missing momentum, which has been previously used by CLEO [6]. 
The requirement on cos 9 t hr reduces QED-related continuum backgrounds by 
50% with a signal inefficiency of 10%. 

The acceptance of the selection requirements as a function of generated q 2 for 




(1) 
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events containing electrons in the momentum interval 2.2 — 2.6 GeV/c from 
B — > X u ev e decay is shown in Figure 1. 
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Fig. 1. Acceptance of the selection requirements as a function of generated q 2 for 
events containing electrons in the region 2.2-2.6 GeV/c from B — ► X u eu e decay. 

3. 3 Backgrounds 

Sources of background for B — > X u ev e include continuum events, hadrons 
misidentified as electrons ("fakes"), decays B — > X c eu e , and various secondary 
decays of 5-mesons. We describe below our evaluation of each as well as 
our procedures for evaluating the associated contributions to the systematic 
uncertainty. 

The continuum contribution is evaluated using the OFF data set. To account 
for the small difference in the momentum spectra due to the beam energy 
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difference (0.6%), the electron momenta in OFF data are scaled by the ratio 
of ON to OFF CM energies. The yields in OFF data are then scaled by the 
factor a = 3.005 ±0.015, determined by the ON to OFF ratio of Bhabha event 
yields in the barrel ECL. As a check of this procedure, we measure the yields 
in the momentum range 2.8 — 3.5 GeV/c, above the kinematic maximum for 
BB events. The resulting signal of 85 ± 93 electron candidates is consistent 
with zero, as expected. We assign systematic uncertainties based on a MC 
study of the detector response to Bhabha events (0.4%), and the discrepancy 
with a as calculated, similarly, with dimuon events (0.2%). 

The remaining contributions to background are from BB events and are esti- 
mated using a large Monte Carlo simulated sample of generic BB events [31] 
that contains roughly three and a half times the number of BB events in the 
ON sample. The MC yield due to fakes from charged pions is corrected for the 
difference in fake efficiency measured by data and MC samples of Kg — > 7r + 7r~ 
decays. The error on the correction (~ 30%) is assigned as the systematic un- 
certainty on the yield from pions. Additional minor contributions from kaons, 
protons and muons to the overall fakes yield are conservatively assigned sys- 
tematic uncertainties of 100%. 

Of the real electrons, those not from primary B — > X c eu e , secondary back- 
grounds, are estimated using the Monte Carlo simulated generic BB event 
sample with electronic branching fractions of D°, D + , J/t/j, ip{2S), D s and 
r assigned according to the current world averages [32]. To avoid any bias 
from a possible mis- modelling of data in MC we use the vetoed J/ip sample 
to measure the normalisation factor for both fake and secondary background 
MC yields. This factor is calculated from a fit of the MC-simulated momentum 
spectrum of vetoed electron candidates from J/ip in 5-meson decays to the 
equivalent spectrum obtained from the data. This sample is statistically inde- 
pendent of the final event sample, and, moreover, may not contain neutrinos 
from primary 5-meson decays, which is the case for events providing fake and 
secondary backgrounds. Contributions from secondary electrons are assigned 
systematic errors based on the electronic branching fraction uncertainties and 
the difference between the MC normalisation calculated as described above 
and that according to the number of T(4S*) events (6%). Overall, the latter 
uncertainty has a less than 0.5% effect on the eventual signal yield. 

The spectrum from B — > X c eu e is modelled using three components: X c = D 
(HQET [33]); D* (HQET [33]); and higher resonance charm meson states D** 
(ISGW2 [28]). To improve the agreement with data, we re- weight the D and 
D* components according to q 2 in order to match spectra generated with world 
average values of form factors [32]. The ratio of D to D* branching fractions is 
fixed according to their world average measurements [32]. The proportion of 
the (D+D*) component with respect to the D** component is determined from 
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a binned maximum likelihood fit [34] of the ON data in the LO region 5 , where 
the B — * X u ev e component is modelled using the inclusive model described 
earlier and fixed such that B(B -> X u ev e ) = (0.25 ± 0.02)% [14]. For the 
D** sub-components, D l and D* we set g(B ^ ( ^+fff^ D ^) = 0.35 ± 0.23, 
which has been determined by averaging maximum and minimum assessments 
of their respective world average branching fractions [32]. Semileptonic spectra 
are also re-weighted to include the effect of QED radiative corrections as calcu- 
lated with the PHOTOS algorithm [35]. It has been observed that the KEKB 
collision energy has variations of (9(1 MeV) over time and that this results in 
a measurable variation of the 5-meson momentum over the running period of 
our ON data sample. As our Monte Carlo generator assumes a fixed energy, 
we apply a shift to the reconstructed momentum in the MC to correct for the 
difference. The correction depends on the beam energy measurement in the 
same run period as our ON data set, which is made using a fully reconstructed 
5-meson decay sample. All spectra for backgrounds other than B — > X c eu e 
are derived from the generic BB MC sample and handled in the same manner 
as for the HI region. The goodness-of-fit as estimated by the x 2 / n df gives 
17.8/13. We use this fit result to not only determine the B — > X c ev e back- 
ground level in the HI region (2.2 — 2.6 GeV/c), but also simultaneously in 
the signal regions defined as 2.3 — 2.6 GeV/c and 2.4 — 2.6 GeV/c. 

The same procedure as described above is repeated for three additional HI 
regions 1.9 - 2.6 GeV/c, 2.0 - 2.6 GeV/c, and 2.1 - 2.6 GeV/c. In each case 
the LO region is adjusted such that its upper bound equals the lower bound of 
the HI region, giving respective LO regions of: 1.5—1.9 GeV/c; 1.5—2.0 GeV/c; 
and 1.5-2.1 GeV/c. The xV ndf for fits in these L0 regions are: 6.8/7; 11.9/9; 
and 13.9/11, respectively. 

The systematic error on our measurement due to the uncertainty in the B — > 
X c ev e shape and relative normalisations is estimated by varying the parame- 
ters fixed in the fit by their individual uncertainties, as described above. The 
maximum deviation observed from either an upward or downward variation 
is assigned as the systematic error. We also calculate uncertainties for cases 
of: no QED radiative correction; an ISGW2 modelled B — > X u ev e spectrum 
shape [28]; and the inclusion of a non-resonant B — > D^neu e (Goity and 
Roberts [36]) decay component in the fit. CLEO included the B — > D^neiy e 
component in their standard fit [6], but in our case, the shape of its mo- 
mentum spectrum bears too close a resemblance to that of the B — > D**eu e 
component. If both D^n and D** components are included in the fit the D** 
component floats to zero. This is clearly contrary to observation, given the 
measured inclusive branching fraction B(B -> D**ev) = (2.70 ± 0.70)% [32]. 



5 The sum of the yields of the components in the fit is constrained to equal the ON 
data yield. 
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The systematic that has the greatest effect on the X c background estimation in 
the HI region is the uncertainty on the D* form factor, which has been obtained 
by varying the form factor slope parameter p 2 within its uncertainty. The net 
systematic uncertainty is calculated as a sum in quadrature of the individual 
systematic uncertainties. Table 1 lists the electron candidate yields in ON 
data, the estimated background contributions and the subsequent extracted 
signal for the six overlapping momentum intervals. 

Table 1 

The B — ► X u ev e endpoint and background yields in six momentum intervals, where 



the first error is statistical and the second is systematic. 



Pcm (GeV/c) 


2.4-2.6 


2.3 - 2.6 


2.2 - 2.6 


Won 


1741 


3534 


8854 


aiVoFF 
Fakes 

N B^J/^e 
NB^X c eu e 


1166 ±59 ±5 
12 ±2 ±4 
40 ± 3 ± 4 
8± 1 ± 1 
4±1±0 


1878 ± 75 ± 8 
43 ± 4 ± 14 
94 ± 5 ± 8 
23 ± 2 ± 2 

345 ±11 ±23 


2743 ±91 ±11 
85 ± 5 ± 28 
191 ± 7± 17 
53 ± 4 ± 4 
3658 ± 36 ± 151 


N B^X u eu e 


512 ± 73 ± 7 


1152 ±97 ±29 


2124 ± 136 ± 155 




PCM (GeV/c) 


2.1 - 2.6 


2.0 - 2.6 


1.9 - 2.6 


Non 


23617 


54566 


104472 


aN FF 
Fakes 

NB^XceUe 


3738 ± 106 ± 15 
93 ± 6 ± 34 
336 ± 9 ± 29 
127±5±10 
15494 ± 73 ± 437 


4900 ± 121 ± 20 
131 ± 7 ± 40 
562 ± 12 ± 49 
263 ± 8 ± 22 
42769 ± 123 ± 970 


6234 ± 137 ± 25 
143 ± 9 ± 47 
880 ± 15 ± 77 
553 ± 11 ±48 
87705 ± 180 ± 1550 


^B^X u eu e 


3830 ± 201 ± 439 


5941 ± 291 ± 972 


8957 ± 395 ± 1553 



4 Extraction of the partial branching fraction 



The inclusive partial branching fraction is found using 
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where N BB = (29.4 ± 0.4) x 10 6 and the overall selection efficiency is cmc- 
The systematic uncertainty on the efficiency includes effects from tracking, 
electron identification, event selection, or model criteria: 

• The uncertainty on the track finding for our electron candidates is studied 
using the MC simulated track embedding method. Care is taken to consider 
all known sources of uncertainty in the MC simulation: magnetic field effects; 
CDC wire hit inefficiency; uncertainties in the material budget of the SVD 
and CDC; and drift time resolution effects in the CDC. The ratio of data 
to MC single track reconstruction efficiency is consistent with unity at the 
1% level. Accordingly, this uncertainty is assigned as the systematic error 
on the efficiency due to tracking; 

• The uncertainty in electron identification (ID) efficiency is measured using 
inclusive J ftp events (The method implemented is similar to that described 
in Ref. [30]). The study involves reconstructing J/ip — > e + e~ decays with 
both tracks satisfying the same track requirements as those of the electron 
candidates considered for this analysis. We find excellent agreement of the 
MC simulation with data at the level of 2% and subsequently assign this as 
the systematic uncertainty on electron identification; 

• The effect of the main event selection criteria, namely those of /C depen- 
dent ^flow an d cos 6 t hr cuts, is assessed in two control samples. We fully 
reconstruct B + — > D (— > K + n~)p + (-^ 7r + 7r°) decays. Here the kaon, disre- 
garding particle identification, is assigned as the electron candidate, whilst 
the pion is regarded as the neutrino. In comparison to B — > X u ev e) the 
mass of the D meson fixes q 2 = m 2 D . The data to MC ratio of the selec- 
tion efficiency is calculated as a function of CM momentum in the range 
1.5 — 2.6 GeV/c, in bins of 0.05 GeV/c; the best fit is achieved with a con- 
stant, which is found to be consistent with unity within 2% uncertainty. We 
also fully reconstruct J/ijj — > e + e~ decays and subtract off backgrounds to 
yield B — > J/ipX decays. We assign the highest momentum electron from 
the J/tp decay to be the electron candidate. The remaining electron is re- 
garded as the neutrino. The mass of the J/ip- meson fixes q 2 = mju. The 
selection efficiency in this sample was measured as described above, and 
the best fit, which was also achieved with a constant, found data and MC 
agree to within 3% uncertainty. Accordingly an overall uncertainty of 4% is 
assigned as the systematic error due to event selection; 

• Model dependence is assessed using the four inclusive samples described 
above. The maximum shift in selection efficiency is assigned as the sys- 
tematic uncertainty due to model dependence, and is dependent upon the 
particular HI region. It varies from 1.7% to 3.4% as the lower momentum 
limit is increased. 

The efficiencies for selecting electrons from B — > X u ev e decays after all se- 
lection criteria have been applied are given in Table 2. Our total efficiency 
decreases as the lower limit of the electron momentum interval increases, an 
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effect due mostly to the momentum dependence of the K, dependent J^flow 
cut. 

Figure 2(a) shows the ON and scaled OFF momentum spectra along with the 
total background. Figure 2(b) shows the ON spectrum after background sub- 
traction and efficiency correction, revealing the contribution of B — > X u ev e . 
The shape prescribed by the inclusive model described earlier, with final state 
radiation, is also shown. The partial branching fractions for each momentum 
interval are given in Table 2; as the lower momentum limit is decreased the 
uncertainty comes to be dominated, as expected, by the uncertainty in the 
B — > X c ev e background subtraction. Our partial branching fraction measure- 
ments are consistent with those of CLEO and have overall reduced uncertain- 
ties [6]. 



5 Extraction of \V u b\ 



5.1 DFN method 



The value of \V u b\ is extracted using the formula given in Ref. [14]: 



\V uh \ = 0.00424^?^|^3^^ x (i.o ± .028 Al , 2 ± 0.039 m J, (3) 

which is an updated version of the expression given in Ref. [32], and includes 
the latest measurements of the heavy-quark expansion parameters [37]. We 
average the current world average neutral and charged 5-meson lifetimes to 
obtain t b = 1.604 ±0.011 ps [32]. To obtain the full inclusive rates for charm- 
less semileptonic _B-meson decay from our partial branching fractions, we must 
determine the spectral fractions /„ and the spectral distortion due to final state 
radiation loss c^rad such that 

B(B -+ X u eu e ) = + w _ (4) 

Ju 



The value of 5rad is calculated from a comparison of MC signal events gener- 
ated with and without the PHOTOS algorithm implemented. It has been the 
convention to assign a 10% systematic uncertainty on the correction based on 
studies that compare the PHOTOS performance with next-to-leading order 
calculations in B — > Deu e decays [38], since the study has yet to be extended 
to B — > X u eu e decays, we assign a larger uncertainty, equivalent to a third of 
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Fig. 2. The electron momentum spectrum in the T(4S I ) rest frame: (a) ON data 
(filled circles), scaled OFF data (open circles), sum of scaled OFF data and esti- 
mated BB backgrounds (histogram), (b) ON data after subtraction of backgrounds 
and correction for efficiency (filled circles) and model spectrum of B — > X u ev e 
decays with final state radiation (histogram, normalised to the data yield in the 
1.9 — 2.6GeV/c momentum range). 



the size of the effect. The correction factors for each HI region are given in 
Table 2. 



The values of /„ for the different momentum intervals are determined with 
the DeFazio-Neubert prescription, using three different forms of the shape 
function with the parameters, A SF and \f F , determined from fits to the Belle 
measured photon energy spectrum in B — > X s 7 decays [21,24]. The resultant 
values of f u are given in Table 2, they range from 3 — 32% as the lower 
momentum limit is decreased. The statistical uncertainty, averaged over each 
shape function form, is determined from the half- difference of maximum and 
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Table 2 

Branching fractions and extraction of |V^&|(DFN method). The reconstruction ef- 
ficiency, eMCj as calculated from Monte Carlo. The partial branching fractions, 
AB u (p), where the errors are from statistics and experimental systematics, respec- 
tively. The lepton momentum spectral fractions, f u , where the first error is the 
combined statistical and systematic uncertainty, and the second error is the theo- 
retical uncertainty in extracting shape function parameters from B — > X s ~f decays 
and applying this knowledge to B — > X u lvi decays. The correction due to the final 
state radiation loss is denoted <5rad- The full branching fractions, B(B — > X u lu{), 
where the first error is due to experimental uncertainty and the second is from f u . 
The extracted values of \ V u b\: the first error is experimental; the second error is from 
f u , combined statistical and systematic; the third error is from f u theory; and the 



last is from the 


application of the V u b formula given 


in Eqn 3. 


PCM (GeV/c) 


€MC(%) 


A^ n (p)(10" 4 ) 


fu 


1.9-2.6 


18.0 ±0.9 


8.47 ± 0.37 ± 1.53 


0.321 ± 0.022 ± 0.041 


2.0 - 2.6 


17.6 ±0.9 


5.74 ±0.28 ±0.98 


246 ± 020 ± 042 


2.1 - 2.6 


17.2 ±0.9 


3.78 ± 0.20 ± 0.48 


0.173 ±0.017 ±0.040 


2.2 - 2.6 


16.6 ±0.9 


2.17 ±0.14 ±0.20 


0.109 ±0.013 ±0.034 


2.3 - 2.6 


16.5 ±0.9 


1.18 ±0.10 ±0.07 


0.058 ±0.010 ±0.025 


2.4-2.6 


16.2 ±1.0 


0.53 ±0.07 ±0.03 


0.025 ± 0.006 ± 0.014 


PCM (GeV/c) 


<5rad 


B(B - X u ev e )(W- 3 ) 


|V; b |(10- 3 )(DFN) 


1.9 - 2.6 


0.06 ± 0.02 


2.80 ±0.52 ±0.41 


5.01 ± 0.47 ± 0.17 ± 0.32 ± 0.24 


2.0 - 2.6 


0.07 ±0.02 


2.49 ± 0.45 ± 0.47 


4.73 ± 0.42 ± 0.19 ± 0.40 ± 0.23 


2.1 - 2.6 


0.07 ±0.02 


2.34 ±0.33 ±0.59 


4.59 ± 0.32 ± 0.22 ± 0.53 ± 0.22 


2.2 - 2.6 


0.09 ± 0.03 


2.16 ±0.25 ±0.73 


4.41 ± 0.25 ± 0.27 ± 0.69 ± 0.21 


2.3 - 2.6 


0.10 ±0.03 


2.22 ±0.24 ±1.02 


4.47 ± 0.24 ± 0.36 ± 0.96 ± 0.22 


2.4-2.6 


0.11 ±0.04 


2.39 ± 0.38 ± 1.46 


4.63 ± 0.37 ± 0.53 ± 1.32 ± 0.22 



minimum /„ found on the A% 2 = 1 contour in (A SF , Xf F ) parameter space. 
The systematic error stems from variation of the scale used to evaluate the 
strong coupling a s (/j) (fi = mb/2,2mb) and differences among shape function 
forms. The theoretical uncertainty is obtained by varying the parameters by 
±20%, reflecting the fact that the procedure is correct only to leading order 
in the HQET expansion that describes the non-perturbative dynamics of B- 
mesons. Our variation is twice that considered by CLEO (±10%). At the time 
of their evaluation little was known about sub-leading and weak annihilation 
effects, they have since been predicted to be large [39,40,41] and are better 
represented by a ±20% variation. The resulting full branching fractions and 
extracted values of |14&| are given in Table 2. All the uncertainties contributing 
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to \V u b\ are summarised in Table 3 for each momentum interval. As expected, 
as the lower momentum cutoff is decreased, the uncertainty from f u that is due 
to theory, decreases, while the main experimental systematic, the estimation 
of B — ■> X c eu e , increases, in line with its background contribution. The best 
overall precision (13%) on \V u b\, based on a sum in quadrature of experimental 
and theoretical uncertainties, is found for the 1.9 — 2.6 GeV/c momentum 
interval with 



V ub \ = (5.01 ± 0.47 ± 0.17 ± 0.32 ± 0.24) x 10" 3 , (5) 



where the first error is from experiment, the second and third are due to 
experiment and theory errors on f u , respectively, and the last is the uncertainty 
in applying the \V u b\ extraction formula. 



Table 3 



Uncertainties contributing to the determination of |V^&|(10 3 ) (DFN method). The 
total error is obtained from a sum in quadrature. 







Momentum Interval (GeV/c) 




Source of Uncertainty 


1.9 - 2.6 


2.0 - 2.6 


2.1 - 2.6 


2.2 - 2.6 


2.3 - 2.6 


2.4-2.6 


Statistical 


0.11 


0.12 


0.12 


0.14 


0.19 


0.33 


B — ► X c eu e background 


0.43 


0.39 


0.26 


0.16 


0.04 


0.00 


Other B background 


0.03 


0.03 


0.03 


0.04 


0.03 


0.03 


Efficiency-detector 


0.12 


0.11 


0.11 


0.10 


0.11 


0.12 


Efficiency- model 


0.04 


0.04 


0.05 


0.05 


0.05 


0.08 


Nbb 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


<5rad 


0.05 


0.05 


0.05 


0.06 


0.07 


0.08 


f u statistical 


0.17 


0.19 


0.22 


0.26 


0.35 


0.51 


f u systematic 


0.04 


0.05 


0.05 


0.06 


0.09 


0.13 


f u theory 


0.32 


0.40 


0.53 


0.69 


0.96 


1.32 


\V ub \ :mf(lGeV), Ai, 2 


0.24 


0.23 


0.22 


0.21 


0.22 


0.22 


Total 


0.64 


0.66 


0.69 


0.81 


1.07 


1.48 
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5.2 BLNP method 

In this prescription |T4&| is obtained directly from the partial branching frac- 
tion A£>, using the formula 

, v , / Agq + ^RAp) 1 

IKbl = V T B r> (6) 

where tb and <5rad are as the same as described previously, and R is the theo- 
retical prediction of the partial rate of B — > X u lv decay in units of | V^ft| 2 ps — 1 
for a given momentum region. The implementation of the BLNP method relies 
on a model for the leading order shape function that is constrained by HQET 
parameters; the mass and average momentum squared of the 6-quark, m&(SF) 
and /z 2 (SF), respectively, as defined in the shape function scheme (SF). They 
are set to m fe (SF) = (4.52±0.07) GeV/c 2 and /4(SF) = (0.27±0.13) GeV 2 /c 2 , 
as derived from the measurement procedure that is described in the introduc- 
tion. 

Table 4 gives the R and \V u b\ values for the overlapping momentum intervals. 
The first error on R is the experimental uncertainty on the leading order 
shape function, which is our own estimation calculated as the half-difference 
of minimum and maximum R values obtained from a set of shape function 
parameters which lie on the A% 2 = 1 contour. The second error on R is a 
theoretical uncertainty stemming from the variation of the matching scales 
fj,i, p, [/,h, sub-leading shape function models and the weak annihilation effect, 
where the latter effect is constant (±1.40) for all momentum intervals [25]. 

Our most precise value, which has an overall uncertainty of 13% as based on 
the sum in quadrature of all the uncertainties, 

\V ub \ = (5.08 ±0.47 ±0.421°;!) x 10" 3 , (7) 

is found for the 1.9 — 2.6 GeV/c momentum interval. When the shape func- 
tion parameters and consequently R are better determined, \V u b\ can be re- 
calculated from the partial branching fraction measurements presented in Ta- 
ble 2. 



6 Summary 

We have measured the inclusive charmless semileptonic S-meson decay branch- 
ing ratio in six overlapping momentum intervals that encompass the endpoint 
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Table 4 

Predicted partial rate R for B — > X u lv\ and extracted value of |V^,| (BLNP method). 
The first error in R is the shape function error stemming from the uncertainty in the 
knowledge of HQET parameters and the second is a theoretical uncertainty stem- 
ming from the variation of the matching scales ji^, p,, fi^, subleading shape function 
models and the weak annihilation effect. The first error in \ V u b\ is the experimental 
error, and the remaining errors are those propagated from R, respectively. 



Pcm (GeV/c) 


i?(|K 6 | 2 ps _1 ) 


\V ub \ (10" 3 )(BLNP) 


1.9-2.6 
2.0-2.6 

2.1 - 2.6 

2.2 - 2.6 

2.3 - 2.6 
2.4-2.6 


21.69 ± 3.62 tlH 
16.05 ± 3.05 li;^ 
10.86 ± 2.51 t\f 7 
6.46 ±1.541^ 
3.15 ± 0.881 lf 4 
1.12 ± 0.391 HI 


5.08 ± 0.47 ± 0.42 
4.87 ±0.43 ±0.461^1; 
4.83 ± 0.33 ± 0.56 Igjjjj 
4.77 ± 0.26 ±0.57±g;|| 
5.07 ± 0.71 ± 0.52 tUl 
5.70 ± 1.00 ±0.67im 



of the electron momentum spectrum. These included a momentum interval 
with a minimum lower momentum cutoff of 1.9 GeV/c, from which the partial 
branching fraction was measured to be AB = (8.47 ± 0.37 ± 1.53) x 10~ 4 . We 
have extracted \V u b\ using both the DFN and BLNP methods, but we adopt 
the results of the latter method since it is more advanced. The most precise 
| V u b | value was extracted from the decay rate in the 1.9 — 2.6 GeV/c momen- 
tum interval and found to be \V ub \ = (5.08±0.47±0.42tg;!|) x 1CT 3 . Owing to 
updated knowledge of background shapes and normalisations, as well as the 
improvement in the theoretical prediction of the decay rates for B — > X u ev e 
and B — > X s ^ decays, the precision of the present measurement is better than 
that of the previous endpoint measurement by CLEO [6]. Although endpoint 
methods have not been preferred for a precison determination of \V u b\ from 
inclusive decays [42,43], the results presented in this letter for the momentum 
interval 1.9 — 2.6 GeV/c are competitive in precision with measurements that 
have utilised the favoured kinematic regions of hadronic mass and dilepton 
mass squared [8,9]. This competitiveness is due to a minimum lower momen- 
tum cutoff of 1.9 GeV/c. Our results also, independent of the extracted value 
of | V u b |, help to bound theoretical uncertainties that in general are encountered 
in all | | extractions from inclusive charmless semileptonic B- meson decays, 
for example, those relating to quark-hadron duality and the weak annihilation 
effect [44]. 

The comparison of our result with other experimental measurements of \ V u b\ [6,8,9] 
must be made on a consistent basis, that is, the extraction of \V u b\ from a par- 
tial branching fraction measurement needs to be performed using a common 
theoretical framework with common inputs. 
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